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Abstract — A novel methodology is introduced here that exploits 2D images of arbitrary elastic body deformation instances, so as to 
quantify mechano-elastic characteristics that are deformation invariant. Determination of such characteristics allows for developing 
methods offering an image of the undeformed body. General assumptions about the mechano-elastic properties of the bodies are 
stated, which lead to two different approaches for obtaining bodies' deformation invariants. One was developed to spot deformed 
body's neutral line and its cross sections, while the other solves deformation PDEs by performing a set of equivalent image operations 
on the deformed body images. Both these processes may furnish a body undeformed version from its deformed image. This was 
confirmed by obtaining the undeformed shape of deformed parasites, cells (protozoa), fibers and human lips. In addition, the method 
has been applied to the important problem of parasite automatic classification from their microscopic images. To achieve this, we first 
apply the previous method to straighten the highly deformed parasites and then we apply a dedicated curve classification method to 
the straightened parasite contours. It is demonstrated that essentially different deformations of the same parasite give rise to practically 
the same undeformed shape, thus confirming the consistency of the introduced methodology. Finally, the developed pattern recognition 
method classifies the unwrapped parasites into 6 families, with an accuracy rate of 97.6 %. 

Index Terms — deformation invariant elastic properties, automatic curve classification, parasite automatic identification, straightening 
deformed objects, image analysis, elastic deformation, pattern classification techniques. 



1 Introduction 

THERE are numerous applications, where bodies suf- 
fer deformation due to elastic forces (stresses). In 
these cases, one frequently encounters two important 
problems: a) to make consistent and reliable estimation 
of the body undeformed shape from images of random 
instances of body deformation and b) to identify the 
deformed body from these images. We would like to 
emphasize that, as a rule, identification of bodies on 
the basis of images of their deformation, is practically 
prohibited by the randomness of the deformation. One 
encounters such problems in various disciplines appli- 
cations, such as automatic identification of highly de- 
formed parasites, cells or large molecules from their im- 
ages obtained via microscope, like the problem tackled 
in (16), in strength of materials, elastography |20|, |9|, in 
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civil engineering in general, etc. 

In the bibliography, there are some approaches for 
generating the phases between an initial and a final 
stage of a body elastic deformation (14) , [19|. In these 
publications, images of both the initial and final stages 
deformation are available. On the contrary, in the present 
paper, images of the undeformed body are not a priori 
available, hence there are no initial phase representa- 
tions. Moreover, methods that do not demand given ini- 
tial body image to perform deformation, use deformable 
lines or surfaces (2), (17), which are not uniquely con- 
nected to deformation invariants of body shape. Thus, 
they cannot be used to create representative undeformed 
body images from its deformed instances. There are, 
also, some approaches dealing with deformable contours 
in the concept of their invariant formulation ||6j, (7). In 
these references, active contour models are formulated so 
as their body shape modeling performance is invariant 
under affine body shape transformations. Namely, in 
[ 7 1 the invariant body shape modeling is formulated by 
finding the eigenvectors of the shape matrix of the body 
contour vectors and in (6) the energy function mini- 
mized along the deformation of the contour model is 
properly selected so that the performed minimization is 
translation and scale invariant. But, in these approaches, 
shape transformation invariance concerns the performed 
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energy minimization along contour-model deformation. 
Hence, in order to use such transformation invariant 
active contour models to obtain undeformed body con- 
tours given their deformed shapes, we should have had 
in advance a reliable estimation of the undeformed body 
shape. Otherwise, the curve to which contour model 
converges and the initial deformed body shape, are by 
no means, necessarily, in a one-to-one correspondence. 
These facts call upon an alternative approach which is, 
for the first time, presented in this paper. 

Hence, in the present paper, we have introduced the 
following approach to tackle the problem of construct- 
ing identifiable images of undeformed body from its 
deformed instances: We estimate elastic deformation 
invariants (curves and sections) of a body suffering an 
equivalent to 2D deformation from an image of it at 
an arbitrary deformation instance. Knowledge of these 
invariants allows for unwrapping /straightening the de- 
formed body image, a fact that in turn permits appli- 
cation of pattern recognition techniques for the body 
automatic classification /identification. We have applied 
the introduced approach to an important and some times 
crucial veterinary problem, namely the automatic identi- 
fication of domestic animal parasites, from their images 
obtained via microscope |18|. These images represent 
the parasites in a state of serious deformation. But the 
applicability of the introduced methodology seems to 
be quite more general, since the assumptions about the 
mechano-elastic properties of the deformed bodies that 
have been stated and adopted, are plausible and quite 
standard. In order to verify this, we have applied the 
introduced methodology to the determination of the 
undeformed state of elastic fibers, protozoa cells and lips 
expressions. 



2 A Brief Description of the Intro- 
duced Approach and the Resulting Ap- 
plication 

One of the major motivations for the present work 
was the fact that, in many important applications, one 
must identify objects on the basis of their deformed 
images. Very frequently, this is a very difficult task; for 
example, in the considered case of identifying parasites 
from images of them in deformed state, obtained via 
microscope, not even the expert parasitologist can iden- 
tify the parasite gender. Thus, the authors developed 
an original methodology to tackle this type of problems 
which includes the following steps: 

1. A set of general assumptions concerning the elastic de- 
formation of the considered body was adopted. A main 
assumption, frequently encountered in elasticity theory, 
is that there are curves and sections whose characteristics 
remain invariant in any deformation stage. 

2. A framework and a set of related equations concerning 
the body deformation process, have been developed, 
on the basis of the aforementioned assumptions. Since, 



many deformation invariant characteristics are associ- 
ated with curves, the equations governing the body de- 
formation have been written in curvilinear coordinates. 

3. The previous analysis gave rise to two different ap- 
proaches for straightening /unwrapping the deformed 
body. In the first approach, a major invariant curve, 
namely the neutral line, is determined together with 
its cross sections. The fact that the size of both the 
neutral line and its cross sections remain the same, leads 
to the construction of the undeformed body shape. In 
the second approach the equations governing the body 
deformation are written in such a manner, so as to corre- 
spond to a sequence of fundamental image morpholog- 
ical operations. Therefore, the image of the undeformed 
body is obtained by application of the inverse procedure 
to the arbitrary deformed body image. 

4. The system that was developed on the basis of the 
above approaches was applied to microscopic images 
of parasites suffering high deformation. The obtained 
images of straightened versions of the same parasite cap- 
tured in different deformation instances, were very close 
to each other, thus confirming the validity and efficiency 
of the introduced approach (Section [7) (Figures [§J [7). 
In addition, the expert parasiologist has been able to 
fully identify the parasite gender from the undeformed 
images generated by the method (e.g. Figures 11 12) . 

5. Finally, the authors proceeded in an automated recog- 
nition of the obtained straightened parasite versions. For 
this reason, they have developed a dedicated partitional 
classification algorithm that employs the Euclidean dis- 
tance of the contours of the unwrapped body images. 
The representative contour curve of each group of par- 
asites has been dynamically estimated via minimization 
of the classification error on a Training Set. The method 
offered an overall success identification rate of more than 
97.6%. 



3 Analysis of the Body Elastic Defor- 
mation 

3.1 Adopted Assumptions Concerning the Consid- 
ered Object's Elastic Properties 

1) All object (e.g. parasite) parts are isotropic, homo- 
geneous and continuous, with a good approxima- 
tion at least. 

2) The static equation of balance holds for the de- 
formed element too (1st order Theory). 

3) There exists a piecewise smooth curve of symmetry 
for the undeformed object, or a curve of minimum 
curvature. We will use for it the term " reference 
curve " . 

4) Object's straight line segments, the cross sections, 
which are initially perpendicular to its symme- 
try curve, remain straight and perpendicular to a 
proper corresponding line after the deformation, 
which is usually called neutral line. The neutral 
line is always coplanar with the reference curve 
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and, hence, the whole deformation can be studied 

in two dimensions. 
5) The generated stresses and displacements along 

the object body are linearly related, namely the 

generalized Hooke's law holds. 
The aforementioned hypotheses allow us to study 
the elastic behavior of the object in two dimensions. 
Consequently, the information extracted from the de- 
formed object images may be sufficient for this study; for 
unwrapping the objects and eventually for automatically 
classifying them. 

3.2 Useful Fundamental Mechano-elastic Quantities 

We now proceed to summarize some notions that are 
commonly used in most approaches of Elasticity Theory. 
These will later be used in the analysis on which the un 
wrapping of deformed bodies is based. In fact, in elastic 

ity theory one defines the stress tensor a = ° xx Gxy 

®yx ®yy 

which expresses forces per unit length in all directions; in 
&xxr c X y, &yx, &yy the first subscript denotes the direction 
to which the stress is vertical, while the second subscript 
denotes the vector component direction. Thus a yx is the 
x-component of the force per unit length exerted on a 
differential element vertical to the y-direction. Similarly 



one defines the strain tensor e 



iyx 



'xy 

z yy. 



which 



expresses the relative elongations along all directions 
and possible torsions. 

The body constitutive equation relates the stress ten- 
sor a with the strain tensor e. These two tensors can 
be related through any functional form. However, in 
many practical circumstances, this functional form can 
be considered to be linear, namely a = Ee, where E is a 
constant matrix (generalized Hooke's law). 

3.3 Body's 2D Deformation Equations 

We assume that the undeformed object has a reference 
curve M : fl(s) = (x^s), y^s)) parameterized via its 
length s. Then we define the unit vector l(s) = yr^jj- 
tangent at an arbitrary point of M at length s and the 
unit vector n(s) normal to l(s). Next, the coordinates of 
the points of object's body will be expressed via their 
distance vector from the reference curve M and the 
length of this curve. Namely for each point with position 
vector xi+yj we express this point via its distance vector 
from M. 

r(s, 5) = fl(s) + (x- x^(s))i + (y - y»(s))j = fl(s) + Sn(s) 

(1) 

Thus any differential displacement in the undeformed 
body is written df (s, S) = (l+5K,(s))dsl(s)+d5n(s), where 
the curvature k(s)1(s) = -^n(s). Now, we can represent 
any differential deformation dui + dwj by using the basis 
(Z(s),n(s)) as a curved vector. 

dh(s) x 



dui + dwj =11 + S UiW l(s) 



ds n 



ds UiW l(s) + d5 u 



h(s) 
(2) 



p = (x=y) 




Fig. 1. Curvilinear coordinates (s,5) of an arbitrary point 
p=(x,y) along a reference line. 



ds UfW = \Jd s u 2 + d s w 2 = dsJ (v T u + (Vw 

(3)_ 

dS u w = ^d 5 u 2 +d 5 w 2 = dS\J(V T u n(s)) 2 + (V T w n(s)) 2 

(4) 

Equivalently, one can define a displacement (d s u + 
d s w)l{s)+d$w, with the functions u and w corresponding 
to u and w via the expressions 



\7u = f( 5 )y (v T u />)) + (v T ^ />)) 

Vw = h(s)yJ(V T u n(s)) 2 + (V T w n(s)) 2 (5) 



Then the differential deformation is written 



(d s u + d 3 w)l(s) + dswh(s) 



|| Vfi|| +1^(5) 



dh(s) 
ds 



It also holds that 

dh(s) d ( Vw 



dsl{s) + \\Vw\\d5h{s) (6) 
H(w) - Iw nn _ t{s) = <s)t{s) (7) 



ds ds \ \\Vw\\J || Vw 

Then u and w are given by the differential equations 

dw 
~d5 



Ww\\ 



dw 
~ds~ 



w„ 



-w 



du 
ds~ 



..Vw|| 
||Vfi|| 



(8) 
(9) 



(10) 



Strains en, s nn , ei n , e n i caused by the deformation 
functions u, w are given by the following expressions: 



£11 



|V«|| 
|Vw| 



1 
1 



= £nl = -. {w -s )— 



(ii) 



According to the Assumption 5, the relation connecting 
stresses and strains along l(s), n(s) is expressed via the 

En E in 



relation a 
matrix. 



Ee; where E 



E n I E n 



a constant 
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Fig. 2. Elementary deformation of a differential element 
npp + n + . This element is placed on the undeformed 
reference line at length s. l(s) is the unit tangent vector 
of the symmetry line at this position and n(s) the unit 
normal one. The differential deformation of this element 
transforms it to nyp /+ n /+ . As mentioned above d s w 
describes differential rotation while d s u and dsw describe 
tangent and normal stretches respectively. 

3.4 Solution of Deformation Functions' PDEs Via 
Morphological Operations 

For an arbitrary point (x, y) we consider the definitions 

A) of its initial distance from a point E of the reference 
line suffering s-deformation obtained via 

{So(x,y)}(s) = nun j yj (x - x^{r)f + (y - ^(r)) 2 | 

(12) 

B) of the correspondence of E with a unique £' of the 
reference line suffering ^-deformation obtained via 

{s (x,y)}(6) = argmin ly/(x - x^(r)f + (y - y^(r)f\ 

(13) 

Then PDE (8) with initial condition {5 (x,y)} (s) is 
equivalent to the action over {5o(x,y)} (s) of a dilation 
filter at scale S with a flat disk kernel (3). Namely, 

w(s,6) = M+[{6 (x,y)}(s)}= sup {{<5 ((a;, y) + e)} (5)} 

Ik1l<<5 

(14) 

Similarly PDE ( [To) with initial condition {s (x,y)} (S) is 
equivalent to the dilation of {s (x,y)} (5) with flat disk 
kernel at scale s. 

u(s,6)=M+[{ So (x,y)}(5)]= sup {{s ((x, y) + e)} (6)} 

m\<s 

(15) 

Finally, equation can be written as In w = 



dn 



(In ||Vw||), with initial conditions {w(x, y)} (0, 5) 



M+[5 (x,y)} and {fey)} (0,J) = ln(||Vt5||) (0,5). In 
Appendix A is shown that the above PDE problem is 
equivalent to constructing the image 



Where the filter a[g(S)](s) is defined by the formula 

a[g(S)} (s) = sup [M s " [<?(£)], inf [M+[g(5)}, ^(s (8))]] 

(17) 

where M \g(S)] — inf {g((x,y) +e)(5)} and the refer- 

M\<s 

ence image {^(x, y)} (s (6)) = -sgn («(s (<S))) oo. The 
— oo normalization is used in order to obtain the identical 
element of inf 



3.5 How the previous results lead to obtaining unde- 
formed body images 

We have reformulated some fundamental relations that 
describe the 2D body deformation, in order to derive 
equations for the deformation process which can be 
subsequently interpreted as a set of morphological op- 
erations acting on a 2D body image. 

Suppose that the unknown image of the undeformed 
2D body is Ijj) evidently, Ijj is unavailable. On the 
contrary, the only available information is an image 
of the deformed body, say Id- We have proved that 
the deformation of the body satisfies equations ([§]>, {9j, 
( [To) . This is equivalent to having proved that Id can 
be obtained from Ijj by application of morphological 
operations on Ijj described in ( [H) , ( [15} , ( (16) . 

Since only Id is available, by applying the inverse 
morphological operation we obtain the undeformed 
image of the body via the process below: 
Initially we create images S = {5 (x,y)} (s ) = 
mm{y/(x - x^{t)) 2 + (y - ^(r)) 2 }, s = 
{so(x, y)} (So) = argmin {(x - x^{r)) 2 + (y - y^(r)) 2 } 

T 

which will play the role of initial conditions for the 
subsequent process: 

A) First we perform dilation of image so until the 
adopted s-scale of deformation and we produce the 
image {s(x,y)}(s) = u[so](s) = M+(s ). Successively, 
we perform dilation of image So until the adopted 
£-scale of deformation and we produce the image 



{s(x,y)} (S) = M+[S ] 



(18) 



B)Next we use this image of (5-deformations to produce 
image 



In 



in I 



(19) 



{w(x,y)} (s,S) = w(0,<5)exp (a 4> (0,8) (s) 



(16) 



on which, subsequently, the following operation is per- 
formed. 

C) We consider that the image {w(x, y)} (s, S) is a 
product of two independent operations {</>(#, 2/)} (s), 

{^(x,?/)} (6). Where {(S(x, 2/)} (5) is the deformation 
that {w(x, y)} (s, S) performs along n(s) and it can be 
immediately obtained from A); while {</>(x, 2/)} (s) is 
the deformation of the body because of the change 
of reference line's curvature. Then, we determine the 
image {(/>(#, 2/)} (s) from the solution of {w(x,y)} (s,S). 
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Namely, {w(x, y)} (s,6) = {s(x, y) } (6) exp (a[g(S)](s)) 
results 

{t(x,y)}( S )=exp(a[g(6)]( S )) (20) 

where g(S) = j(/>o(#,2/)} (8) has been obtained by |l9| . 

Next, on the basis on Section [3] analysis, two ap- 
proaches are analytically described in order to construct 
unwrapped body versions from its deformed images. 

4 The First Approach - Determination 
Of Deformed Body's Neutral Line 

In the present section we will introduce a novel method 
for the determination of the neutral line in the deformed 
body; we once more emphasize that along the neutral 
line no stress is exerted and that this curve corresponds 
to the reference line of the undeformed body, as we will 
demonstrate below. 



4.1 Confirmation that the neutral line passes from 
the middle of its cross sections 

Consider a section A normal to the reference line in the 
undeformed object's state. Then, the vertical shear force 

V and force N normal to A and along the reference line 
are given by the expressions: V = § A cri n dA and V = 
§ A cr u dA 

We assume that the body is in an equilibrium posi- 
tion in each image of its deformation instances. Then, 
equilibrium along reference line implies that N = 0, 

V = 0. Since we have adopted the assumption that the 
generalized Hooke's law holds, it follows that 



N 



= En(\\Vu\\-l)A- 
= Eu(\\Vu\\-l)A 



'in ^nn 



Nw\\ 



(w — So)dA 



(21) 



Combining pi) and N = 0, we deduce that 

{s(x,y)}(6) = s + s (22) 

indicating that s-deformation is only offsetting the length 
of the initial reference line and hence that it can be 
ignored. Now, reference line's stresses in the direction 
tangent to it, an, are given by 



J A 2 ||Vw|| J A 



(23) 

(w - S ) 
(24) 



= E ln £(\\Vw\\ -i)dA 
Combining l |23) and V = 0, we deduce that 

{5(x,y)}(5) = 8o + 6 (25) 

which similarly results that ^-deformation is only offset- 
ting the distances from the undeformed reference line 



and hence that it can be ignored. Now, reference line's 
stresses in the directions normal to it, a\ nr are given via 



(26) 



Combining ( [26] ) and ( [23} with zero offsetting of 5o(x,y) 
we obtain a a = -\ {<j>(x, (s)£ ^gp a Wl 

To study the stress that reference line suffers we 
let £ 0. Then (19 1 implies {<[>o(x,y)\ (5) 

} im J n (mi S o(x,y)) = ]xm \<t>(x,y) 

1 and finally {20} results in lirn^ {w(x, y)} (s, S) 



(s) 



lim So \(f)(x,y)> (s) = 0. Therefore, the curve to which 

<5o— >-0 I J 

symmetry line is transformed suffers no stress; the curve 
formed by all these unstressed points, usually called the 
neutral line, has the following properties: 

A) It is the curve to which the reference line is trans- 
formed due to the elastic deformation process. 

B) No stress is exerted along it. 

C) As a consequence, the neutral line and the body 
reference line are of the same length. 

D) The neutral line passes from the middle point of 
each cross section in the 2D image representation of the 
deformed body. 

E) The undeformed cross sections initially perpendicular 
to the reference line, remain perpendicular to the neutral 
line even after the body deformation. 



4.2 An introduced important property of the de- 
formed object contour tangents 

In this paragraph we will state a new lemma that will 
play an important role for the determination of the 
neutral line and the unwrapping of the body. The proof 
of this lemma is given in Appendix B. 

Lemmma 

Let be the symmetry line of the unwrapped body 
and AD an arbitrary cross section, intersecting at 
point M, where part EiM of the symmetry line has 
length s. Then, due to the deformation, AD moves to 
a section A'D' , perpendicular to the neutral line at point 
M' . Let moreover U and L be the upper and lower 
boundary curves of the deformed body, respectively 
and M'N' be its neutral line differential element. If 
T% is the tangent vector of U' at A! , ffj, the tangent 
vector of V at D' and L M > is the tangent vector of the 
deformed symmetry line at point M'; then it holds that 
Z 



stormed symmetry line at point M ; men it hole 



This lemma proves to be fundamental in defining and 
obtaining the cross section's coordinates, as well as those 
of the neutral line. Knowledge of the coordinates of the 
neutral line allows us to "unwrap" the outline of the 
deformed body, as is shown in Sects. 4.4 4.5 
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Fig. 3. Interpretation of the Lemma. The arbitrary cross section AD in the undeformed parasite body moves to A'D' 
via the deformation process. U' denotes the body upper contour and V the lower contour. is the tangent vector of 
U' at A\ the tangent vector of V at D' and L M ' is the tangent vector of the deformed symmetry line at point M' 
of the deformed neutral line. We determine body differential element's cross section via the demand that £ + = tt. 
Notice that in a neighborhood of D' there is no other point satisfying £ + = it since remains fixed and £ changes. 



4.3 Object Contour Extraction And Its Polynomial 
Approximation 

As the introduced methodology makes use of the border 
line of each parasite, it is necessary to obtain well- 
defined boundary lines of all instances of the examined 
parasites. In order to achieve this, the following method 
is used: 

First, we have applied various image segmentation 
methods ( (T), (8)) in order to obtain quite clear-cut and 
accurate region borders of each parasite. A rather simple 
method that seems to work well, for the considered 
images, is the one that uses each parasite's pixel intensity 
histogram and the lower turning point of it. All pixels 
with intensity lower than this turning point may be 
considered to belong to the parasite body. This method 
may generate various artifacts that may be removed by 
application of proper morphological filters, see e.g. (12) . 

As it will become evident from the subsequent analy- 
sis, in order that the introduced methodology is applied, 
each contour line must have the following properties: a) 
each pixel must have exactly two neighboring pixels, b) 
no isolated pixels or groups of pixels are allowed and c) 
three pixels must not form a compact right (90°) angle. 
Since no edge detection algorithm can generate the 
parasite contour in this form, suitable software has been 
developed to achieve that. The next step is to determine 



if there are specific mathematical curves that optimally 
fit the object body, e.g. a parasite, contour in the obtained 
images. There are various techniques for achieving this 
goal (e.g. (2), (5J, (l3j). For the present application the 
following method proved quiet satisfactory: 

The curve parameter is chosen to be its contour length 
s, calculated via the distance of the successive pixels 
that form it. Subsequently, we approximate the variables 
x and y of the body contour by polynomials of the 
appropriate degree, estimated by means of the body 
contour curvature: 

N N 

x(s) = *ns n , y(s) = h ns n (27) 

n=0 n=0 

Suppose now that one wants to test if the upper body 
contour corresponds to the aforementioned curve and 
at the same time to compute the parameters of this 
optimal curve: Let rf, i = 1,2, ...,N B be the centers 
of the pixels forming the upper contour and let II = 
[a n , ...,ao,b n , ...,6o]- Hence the parametric vector equa- 
tion of the prototype curve is f M (s|II) = x{s)i + y(s)j. 

Next, we compute the optimal set of parameters 11° 
and the corresponding sequence of values of the inde- 
pendent variable s it i = 1,2, ...,7V B , so that r M (s;|n°) 
best fits rf according to the chosen norm E 2 = 

Z^i=i II' i ' i ||- 
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4.4 The Developed Algorithm for Determining the 
Deformed Body's Neutral Line 



In this section, we will analytically describe the method- 
ology we have introduced and applied for determining 
the exact position of the neutral line in the deformed 
body image, as well as the positions of the cross sec- 
tions that, by hypothesis, always remain undeformed 
and normal to the neutral line. The application of this 
methodology has been made on the available parasite 
images and comprises the following steps: 
Step 1 - We, first, extract the parasite contour (Sect. 4.3 >. 




Next, we spot the head and the tail of the parasiteas 
follows: First, to spot the tail, for each contour pixel p 
we consider the sets of pixels Pl that lie on its left and 
Pr that lie on its right. We approximate both Pl and Pr 
with line segments in the Least Squares sense. We let the 
tail T be the pixel where these two line segments form 
the most acute angle. 

Second, we spot the parasite "head" H: We move 
away from the tail and we locally approximate the 
contour by polynomials of fifth degree, of which we 
compute the curvature. We let the "head" be the point 
of maximum curvature which also lies between 0.4 and 
0.6 of the whole contour length. 

Step 2 - We divide the whole contour into two parts I 
and II (arbitrarily upper and lower), that both end at 
the parasite "head" and "tail". Then we approximate 
both parts with polynomials of type ( [27) . All performed 
experiments indicate that this approximation is excellent. 
We form a dense sequence M 1 - 1 , j = 1...N 11 of points ^^2,j = 



Fig. 4. Determination of the cross sections in a highly 
deformed parasite image. A first example. 




Fig. 5. Another example of cross sections determination 



7r| which occurs at point 



Ni 1 , 



belonging to the polynomial curve best fitting part II 
and a less dense sequence M 7 , i = 1...N 1 on the curve 
fitting part I; let r\ and r- 1 be the unit tangent vectors 
to these model curves at each M/ and M 77 respectively. 
Step 3 - Subsequently, we spot parasite's neutral line 
by applying the aforementioned lemma as follows: We 
move away from tail T along part I and we connect 
M{ with each point of set M 77 , j = 1...K, where K 
is a predefined number of pixels, say 5% of the whole 
contour length. We form vectors f\j = M 77 M/. We 
keep only those vectors fij that lie entirely within the 
parasite body and for these we compute the angles <p[j, 
(p[j formed by each vector fij and the tangent vectors 
ff and r(j respectively. Then, we define the sequence 
Acpij = \ip[j + ip[j — 7r| and we let TV 77 be that point 
where the minimum value of the sequence Acpij occurs, 
say the <ii-th of the sequence M 1 - 1 ) we consequently 

define M/iV/ 7 to be a cross section of the parasite that 
remains undeformed and normal to the neutral line. 

Next we compute the second cross section as follows: 
We move away from the tail vertex T, and M{, at 
M\ and once more, we define the set of points Mj^. 
j = 1...K. Proceeding as before we define vectors r^j = 
M^ +j MJ{, j = 1,...,K. We compute the corresponding 
angles ip^j between f 2 j and ff, as well as cp^j between 
7*2, j and ^dl+j- We spot the minimum of the sequence 



say the <i 2 -th point of sequence M 77 . We let M^N^ 1 be 
the second cross section that remains undeformed and 
normal to the neutral line. 

Finally we proceed in obtaining all cross sections 
M/iV/ 7 passing from M 7 , i = 1...N 1 by the same method 
(Figures [3J 5). The middle points of these cross sections 
belong to the neutral line and the unit vector normal to 
these sections is tangent to the neutral line. 

4.5 Unwrapping The Deformed Object 



We have shown in Section 4.1 that under the adopted 
assumptions the neutral line undertakes no stress and 
it is found in the middle of the corresponding cross 
sections. Consequently, we define the neutral line of 
the deformed parasite to be the locus of the middle 
points K v of the cross sections M^N^ 1 as they are 
determined above. Clearly, the length of the neutral line 
remains unchanged during the various phases of the 
body deformation and it coincides with the length of 
the undeformed body's symmetry axis. 

Therefore, in order to straighten the body and find 
its undeformed shape, we proceed as follows: First, we 
compute the distance of all successive middle points Ki 
and Ki+i. Then, along the x-axis we form a sequence 
of points A{ of equal number with Ki as follows: Ai 
is placed in the axis origin; A2 in the positive x-axis, so 
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Fig. 6. Contours of the unwrapped versions obtained 
from six different instances of larva deformations of the 
same individual 
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Fig. 7. Contours of the unwrapped versions obtained from 
three different instances of another individual. 



that the distance between points Ai and A 2 is equal to the 
distance between middle points K\ and K2 . We continue 
this process so that A$ and A i+i are equidistant with 
KiKi+i, until all middle points are exhausted. Moving 
in a direction perpendicular to the x-axis at each point 
A^ we choose two points: Ai with ^/-coordinate equal to 
Xi/2, namely half the length of the cross section M/iV/ 7 
and point B { with ^/-coordinate —Xi/2. We set points A { 
to form the one part of the parasite, while points Bi form 
the other parasite part (Figures [§J [7). 



5 The Second Approach - Determining 
The Undeformed Body Shape Via Filter- 
ing Operations On Its Deformed Image 

The method introduced in this section utilizes the so- 
lution of 2D body deformation PDEs as it was formu- 
lated in Section |3.5| in order to construct the inverse 
deformation for the considered deformed body image. 
Also, exploitation of Assumption 2, in Section |4.1[ has 
simplified this solution, since there, it has been shown 
that solution's Step A can be ignored. So, the sequence 
of image operations that construct the unwrapped body 
version is described in the following subsection. 



Fig. 8. Image, where the parasite body pixels have 
intensity S (x,y) 



5.1 Description Of The Proposed Method. Applica- 
tion To The Unwrapping Of Deformed Parasites. 

Step 1: Initially, the vector parametric equation of the de- 
formed body contour is approximated via polynomials 
of type (27) in the least squares sense, thus obtaining 

f c (s) = (x c (s),y c (s)). Using this vector parametric equa- 
tion, we compute the unit directional vectors along r c (s), 
feOO^OO) - i=M*L*ct*ll , as well as the 



l(s) 



^x c {s) 2 +y c (s) 2 

curvature c{s) = 



, h(s) 



^x c (s) 2 +y c {s) 2 ' 



(* 2 c +y 2 c) 3/2 ' 

Step 2: Afterwards, at any pixel point (x, y), inside the 
deformed body, we attribute the values the follow- 
ing two images, s$(x,y), 8$(x,y), have at this point 

S (x,y) = mm{y/(x - x c (s)) 2 + (y - y c (s)) 2 j, s (x,y) = 

argmin — %c(s)) 2 + (y — 2/c(s)) 2 j. These two im- 

ages, so(x,y), 5 (x,y), play the role of initial conditions 
for the PDEs d8b, JTol, J91 which describe the deforma- 



©/ 0/ © which d( 
tion functional {s(x, y) (5), {s{x, 
respectively. 

Step 3: We have proved in Section 
of these deformation differential 



y)}(s), [j>{x,y)}(s) 



3.4 



alent to applying dilations (10) and ^ 

a[g,K](s 



that the solution 
equations is equiv- 

to s (x,y), 

5o(x,y) respectively and the filter a[g,K](s) to the di- 
lated So(x,y) (see Figure [8 > as formula (|9) indicates. As 
it was shown in Section 4~T] the crucial operation is 
curvature deformation j ?/)j (s) given by formula 

|20l. Namely {<f>(x,y)} (s) = exp(a\g,K 4> ](s)), where 

fy(z,y) = -sgn(c(s (x,2/)))oo, g(x,y) = ln<t> (x,y) and 
cf)o(x,y) = \\\75o(x,y)\\. This image (Fig. |9J) represents all 
curvature deformations applied on the body, through 
scale s. Since we want to straighten parasite bodies 
we demand zero curvature for the reference line. Thus, 
for each pixel (x,y) in the deformed parasite body, 
we determine the proper scale cr(x,y) that minimizes 
{<j>{x,y)}(<r), namely a(x,y) = arg min {{<j>(x, y)} (s)} 

s 

until a curve of a threshold curvature deformation is 
spotted (see Figures |9| (To). Consider that this curve con- 
sists of points y^). Finally, if the given image of the 
deformed body is f(x, y), then the above procedure gen- 
erates an image f T (x,y) = f(x,y), x = M+^ y) [s {x,y)), 
y = S (x,y) - 5 (x fl ,y fl ) (Figures [IT} [l2j. 



Fig. 9. Image of parasite body curvature deformation 
minimized at a scale a(x, y) 




(c) 

Fig. 13. Segmented microscopic images of three de- 
formed elastic fibers of the same type and structure. 



(a) 



(b) 



Fig. 10. Parasite body reference curve of approximately 
zero curvature deformation estimated after inversion. 




(c) 

Fig. 14. Straightened versions of the fibers of Fig. [13] 
presented in the same order with the segmented images. 



Fig. 1 1 . Generating the unwrapped version of the parasite 
of Fig.g 



The case of cell undeformed shape restoration can be 




Fig. 1 2. Generating the unwrapped version of the parasite 
of Fig.|5| 



5.2 Application of the method to the determination 
of the undeformed shape of other bodies 

We have applied the method described above to three 
different cases of body deformation, namely the defor- 
mation of cells (protozoa), in human lips deformation 
and in elastic fibers deformation. 

In fact, the case of elastic fiber deformation (fig. [13} 
can be tackled by immediate application of the method 
described in Sect. |5.1| above, since the fibers have a 
symmetry axis in their undeformed state. We notice that 
both methods described in Sections [2J |5.1| can be applied 
for the determination of elastic fibers undeformed shape 



treated by means of the method described in Sect. 5.1 
since in the undeformed state the cells bear a reference 
curve, which is not necessarily a symmetry axis. In other 
words, Assumptions 1,2,4,5 stated in Section 3.1 hold 
intact both in this and in the elastic fibers case; concern- 
ing Assumption 3, we note that in the case of the fibers 
the reference curve is also a zero curvature symmetry 
axis, while in the cells case this condition is relaxed. We 
note that for the restoration of the undeformed shape 
of the cells considered here, we have employed the 
exact method described in Sect. 5.1 above, using elliptical 
coordinates for the description of the undeformed body 
points (Fig. (16) . 

The case of lips deformation is, in a sense, closer to 
the cells case. In fact, it is quite logical to adopt the 
assumption that there is a reference line of minimum 
curvature in the undeformed lips state and as a conse- 
quence we may apply the method of Sect. 5.1 above. 
We, once more, note that all points of the undeformed 
body are expressed in elliptical coordinates with respect 
to this reference line. The validity of the aforementioned 
assumptions is supported by the following figures. 



10 








(c) (d) 
Fig. 16. Intermediate steps of the deformation inversion 



process of Sect.[57T|for the individual of Fig. |1 5e 

a) lnitialization of the filtering process; image, where the 
cell body pixels have intensity 5 (x,y). 

b) lmage of cell body curvature deformation minimized at 
a scale a(x, y). 

c) Cell body reference curve of minimum curvature defor- 
mation estimated after inversion. 

d) Cell body points registration along reference curve via 
dilation of s (x, y) at a scale of a(x, y) and transformation 
of the resulting registration to elliptic coordinates. 



(a) 





(c) (d) 

Fig. 17. Undeformed versions of the protozoa of Fig. [15] 
presented in the same order with the segmented images. 




(a) 


(b) 














(c) 


(d) 



Fig. 18. Cropped and segmented human lips images of 
two different expressions of the same person. 



6 Automatic Classification of 
Deformed Objects - Parasites 

6.1 Parasite Samples Acquisition 

Diagnosis of the presence of parasitic infections in 
domestic animals, requires a coprological examination, 
where eggs of the parasites strongyle family need to 
be identified. Automatic identification of these parasites 
essentially reduces economic losses due to massive infec- 
tion of domestic animals and it helps rapid restriction 
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(c) (d) 

Fig. 19. Intermediate steps of the deformation inversion 
process o f Sect.[5j]applied to the lips expression of Figs. 
fTSc|[TSa| 

a) lnitialization of the filtering process; image, where the 
lips pixels have intensity 5 (x,y). 

b) lmage of the lips curvature deformation minimized at a 
scale a(x, y). 

c) Lips reference curve of minimum curvature deformation 
estimated after inversion. 

d) Lips points registration along reference curve via dila- 
tion of s (x,y) at a scale of a(x,y) and transformation of 
the resulting registration to elliptic coordinates. 




(a) (b) 



Fig. 20. Versions of the undeformed shape for the two 
different lips images of Fig. [18] It is evident that the in- 
troduced method, applied to the two different expressions 
furnished very similar undeformed states of the lips 



of the disease. The veterinary expert obtained parasite 
images from cultured faecal samples collected from natu- 
rally infected grazing sheep. According to the expert par- 
asitologist, three to six images of each individual motile 
live larva, at lOx objective magnifications, were recorded 
by using a suitable charged coupled device camera 
mounted on the light microscope. A total of 317 images 
of 82 individual larvae in various shapes belonging to 
the genera Trichostrongylus(75 instances), Oesophagos- 
tomum (50 instances), Cooperia(35 instances), Osterta- 
gia(33 instances), Haemonchus(52 instances), Teladorsa- 
sia(72 instances) were recorded. Each image was labeled 
with the genus name of the corresponding larva, its 
magnification and an identification number. 

6.2 The developed method for automatic classifica- 
tion of parasite contours 

We, once more, emphasize that the unwrapping process 
does not take into consideration at all that the deformed 
objects are parasites. Similarly, the automatic classifica- 
tion method presented below is very well applicable 
to any straightened object contours. In both cases the 
only restriction is the validity of the general assumptions 
made in Sects. H] and |U We divide the available set of 
unwrapped parasite curves into a Training and a Test Set 
almost equal in number, by means of a random number 
generator. We also enumerate the individuals of the 
Training Set randomly The expert, e.g. the parasitologist, 
classified all items of the Training Set into a number 
of groups based on the information the unwrapped 
parasites furnished. In order to test the efficiency of the 
classification method with maximum possible reliability, 
we have randomly generated a large number of pairs of 
Training and Test Set, say 200, and we have applied the 
classification method described below for each such pair 
separated. 

6.2. 1 The developed method for automatic classification 
of parasite contours 

First, we arrange all straightened contours of a group 
as follows: We consider the unwrapped interpolated 
curve of the arbitrary i-th element of group k, gener- 
ated by means of the method described in Section E] 
40) = (4( S )>2/U S )) - Ms) = li(s)h k (s), where p k {s) 
the symmetry line anf hk(s) the unit vector normal to its 
tangents h(s). 

We choose the contour of maximum length, say of 
group k with object index M, and we fit the unwrapped 
contours of the other samples of the group to it, via the 
following process: 

Let the contour of maximum length be c^f(s) with 
the independent parameter of the length of the corre- 
sponding symmetry line s £ [0,Ljf]. Then, the other 
objects' contours move along group's symmetry line so 
that their Euclidean distance from cjf(s) is minimal. 
Namely, for the i-th sample of group k, c^(s), s £ [0, L\], 
we choose the translation d\, < d\ < L^f — L\ 
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along group's symmetry line that minimizes the Eu- 

T 1 ii • 1 1 9 



clidean distance XL=o ||^l( 5 ) — ^k( s + 



lently maximizes the quantity Xlf=o 
Then we define the curve that represents /c-group's object 
contours to be the mean curve of group's elements: 

400 = £E£i w 4(*-4)> s £ [o,if], the 

number of /c-group elements with translation d\ > s, 
Nk the number of the elements of /c-group. 

6.2.2 Tuning of the group representative curves of the 
Training Set 

We test the representability of the mean curves generated 
above, for the elements of the Training Set. In fact, we 
first attribute each element of the Training Set, whose 
curve is Cj(s), to the group with mean curve Ck*(s), de- 
termined by the requirement that Euclidean distance of 
the two curves, computed as in Sect. |6.2.1 is minimum. 

When this attribution is over, we check its validity, 
taking into account the parasitologist's classification for 
the Training Set. If and only if a misclassification has 
occurred, we re-estimate the corresponding groups' rep- 
resentative curves as follows: 

Let cj,(s) be the representative curve if the group 
to which (P(s) has been erroneously attributed, while 
let (? r (s) be the representative curve of the group to 
which (?(s) should have been attributed. Then , we first 
optimally match <?(s) to (? w {s) as in Sect. 6.2.1 and after 
that we re-estimate (? w (s) via the formula 

<{s) <- ^(«) + se[o,i&] 

( K 

fa = exp — J2\\c? w (s)-<?(s + di)\ 
Ni <- Ni-Pi (28) 



or equiva- 

cT(* + 4). 



Similarly, we let <?(s) optimally match to (? r (s) and we 
re-estimate c? r (s) by means of the formula 

( -1 Li 

Pi = exp + 

N 3 r <- N 3 r -pl (29) 

In both these re-estimation formulae quantities f3 3 w and 
P 3 r are chosen so as to "punish" large distances between 
element and group curves and at the same time to 
support small distances. 

We apply this process to all erroneously attributed ele- 
ments of the Training Set, thus obtaining a new ensemble 
Eff R of representative curves for each group. 

We repeat the above process by attributing all curves 
(?(s) of the Training Set elements to the representative 
curves of group and we check if the resulting clas- 
sification has a success rate above a certain satisfactory 
threshold, say 98%. If it does we stop the process and 
we consider the elements of E^ R as actual representative 
curves of the existing curves and genera in the consid- 
ered Training Set. Otherwise we repeat the previously 
described re-estimation for the curves of E^ R , thus 
creating a new ensemble E^ R of representative curves 
for each group and so on, until the desired high success 
ratio is achieved. In the end of this process, an ensemble 
E G R of curves c k *(s) is obtained which represent the 
groups of the training set optimally. 

6.2.3 Classification of the elements of the Test Set 

After generating E GR , we proceed in the classification 
of the Test Set parasite images with a method analogous 
to the one described in the previous paragraph. More 
specifically, in connection with each available parasite 
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(a) Plot of all ensembles E GR of groups k=1 and 
k=6 representative curves 

Fig. 22. Training process for the groups k=1 and k=6 representative curves 



(b) The final representative curves of these 
groups, in E GR , that offer success rate over 98%. 



image of the Test Set, we unwrap the deformed par- 
asite by the method introduced in Sect. [2| forming a 
set of straightened parasite interpolated contours, say 
C TS . Next, we compute the distance of each curve &>(s) 
belonging in C TS , from each group representative curve 
of ensemble E GR via the formula 

A(?,4) = ^|^'( 5 )-4( 5 + 4)| 2 , c k (s)eE GR (30) 

s=0 

where d° k is the translation along Ck(s) that centers (?(s) 
and Cfe(s) optimally in the least squares sense. Then, 
we attribute (P(s) and the corresponding parasite to the 
group of minimum distance; namely j-th element of Test 

argmin{A(#',c fe )}- 

k 



Set belongs to the ra-th group if m 



7 Evaluation 
Methodology 



of the Introduced 



7.1 Evaluation of the method of unwrapping the de- 
formed bodies 

If the assumptions made in this paper and the intro- 
duced methodology are correct, one expects that dif- 
ferent instances of a specific body deformation will 
generate very similar undeformed body shapes. In par- 
ticular, one expects that different images of the deformed 
parasite must offer quiet close unwrapped versions after 
application of the methodology introduced in Sections [I] 
and [5] In fact, Figures [6j [7 demonstrate that this is indeed 
the case: The undeformed parasite borders have a differ- 
ence that might be considered negligible in respect to the 
parasite dimensions. We have employed five different 
measures to describe the differences between the shapes 
of the unwrapped parasite that resulted from different 

phases. These measures are: = 1 j P 100%, where q 
is the length of the unwrapped parasite obtained from 
the z-th wrapped larva phase and l p its mean value, 



&2,i = L gp — 100%, where Ef is the area of the un- 
wrapped parasite obtained from the i-th wrapped larva 

mean( yi , j-yj) 

phase and E p its mean value, — — me an(y ) — 100%/ 

3 

where y^j is the width of the unwrapped parasite that 
corresponds to the i-th phase at jDoint xy, hence, we 
define yj = mean^j), a 4 ^ = nt n " P n 100%, where Ilf is 

the perimeter of the unwrapped parasite obtained from 
the i-th wrapped larva phase and fl p its mean value, 

a 5,i = Ci &p 100%/ where Cf is the maximum cross 
section diameter of the unwrapped parasite obtained 
from the i-th wrapped larva phase and C p its mean 
value. 

The mean value and standard deviation of quantities 
&i,i/ &2,i/ o>3,i, o>4,i, are shown in Table [T] 

Very similar undeformed shapes have been obtained, 
when the introduced methods have been applied to 
the images of deformed elastic fibers, protozoa and 
lips expressions. This is shown in Figs. [l4j [Tz] and 
[20| Evidently, the introduced method for deducing the 
undeformed body version from its deformation images, 
seems consistent, reliable and robust. 



7.2 Evaluation of the method of parasites automatic 
classification 



As described in Section |6.1| the entire set of data con- 
sisted of 317 parasite images that have been shot at 
an arbitrary deformation instance. We have randomly 
divided this set into 200 pairs of Training and Test 
Sets almost equal in number. For each Training Set we 
have applied the procedure described in Section |ij thus 
classifying all its members into groups, which in turn 
belong to one of the 6 genera. 

Next, we have considered all samples of the corre- 
sponding Test Set and we have automatically classified 
them to a corresponding group and gender. We have 
repeated the aforementioned process for all 200 pairs of 
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TABLE 1 

Parameters to Measure Consistency of the Unwrapped Parasite Contours 





Evaluation of the discrepancy 


Evaluation of the discrepancy 


Average discrepancy of all 




between the 


unwrapped 


between the 


unwrapped 


experimental results 




versions of Figure 6 


versions of Figure 7 








Mean 


Standard 


Mean 


Standard 


Mean 


Standard 




value % Deviation % 


value % Deviation % 


value % Deviation % 


a i (Length) 


0.82 


0.29 


1.04 


0.73 


0.90 


0.61 


(i2 (Area) 


2.30 


1.18 


1.41 


0.62 


1.67 


1.04 


a 3 (Width) 


1.87 


1.06 


1.32 


0.60 


1.69 


0.91 


<24 (Perimeter) 


0.84 


0.34 


1.07 


0.74 


0.98 


0.61 


a 5 (Max cross section) 


2.63 


1.68 


2.62 


1.63 


2.68 


1.70 



TABLE 2 

Parasite Classification Results for All 200 Randomly 
Generated Test Sets. 



Family 


Number 


Percentage of 
Test Sets with 
Classifications 
of erroneous 
identifications 


Percentage of 
Test Sets with 
Classifications 
of 1 erroneous 
identifications 


l)Cooperia 


17 


53% 


41% 


2)Oesopha- 
gostomum 


24 


30% 


50% 


3)Ostertagia 


16 


100% 


0% 


4)Tricho- 
strongylus 


31 


32% 


48% 


5) Haemonchus 


21 


29% 


71% 


6) Teladorsasia 


29 


52% 


44% 



Training and Test Sets and we have kept record of the 
correct and erroneous test classification of the Test Group 
in each case. The related results are summarized in Table 

m 

7.3 Some aspects concerning the efficiency of the 
introduced methods 

The efficiency of the methods presented here, developed 
for constructing undeformed body shapes and images, 
relies also in the accuracy of segmentation and body 
extraction results. Once the body shape has been ex- 
tracted from the image of its deformation instance with 
noisy but non-conceptual distortion (zero mean value 
noise), both methods developed for constructing the 
unwrapped body version can circumvent this inaccuracy. 
Namely, for both methods described in Sect. 4.4 4.5 and 
in Sect. |5j local distortion of the body contour does not 
affect the polynomial form optimally fitted to it, since we 
expect zero mean value for this local distortion and the 
polynomial fits the contour points in the least squares 
sense. The method employed here for segmentation of 
parasite images makes use of the color intensity local 
homogeneity under the assumption that this homogene- 
ity is represented by a normal distribution of color 



intensity and hence local inaccuracies have zero mean 
value. So, provided that this method offered reliable 
extraction of body shape, local inaccuracies do not affect 
the introduced unwrapping methods' performance. In 
any case, if in another application this segmentation 
method couldn't offer reliable body shapes, different seg- 
mentation approaches (active contours (6), watershed- 
based algorithms |11|, etc.) could be employed. 

Another aspect about the efficiency of the methods 
introduced here, concerns their applicability to the re- 
construction of body undeformed shape in other entities 
than parasites. In the sense that the considered bodies 
suffer an equal to 2D deformation, analysis of Sect. |3.3| is 
valid. But the reliable estimation of the body unwrapped 
version depends on the demand that body deformation 
is invertible and on the reliability of the estimation of 
body undeformed reference curve. Namely, if the as- 
sumptions described in Sect. |3.1| hold for the considered 
body then the performed experiments show that both 
methods described in Sect. |4.4| |4.5| and in Sect. |5.1| are 
directly applicable provided that there is a symmetry 
line in the undeformed body state. This is the case for 
the elastic deformation of various objects /entities, such 
as fibers, cables, nails, sheets, but also serpents and eels, 
etc. On the other hand, there is a variety of entities that 
may suffer elastic deformation of the form described 
in Sect. 13.11 which in their undeformed state have a 
piecewise smooth curve of minimum curvature, i.e. the 
reference curve. Such bodies are cells, viruses, the human 
lips, etc. One may, also, tackle the problem of retrieving 
the undeformed shape of such a body by applying the 
method described in Sect. 5.1 Characteristic examples 
of the aforementioned cases, additional to the parasite 
one, are given in Sect. 5.2 supporting the applicability 
of the introduced methods. We would like to emphasize 
that the introduced methods are valid both for rectilinear 
and curvilinear coordinates describing the body points. 

Concerning the developed curve classification method, 
the efficiency of its classification performance depends 
on how representative are the mean curves of each 
group for the curves that belong to it. In these cases of 
bodies, such as parasites, cells, fibers, sheets, industrial 
line's objects, etc., undeformed shapes of objects of the 
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same kind are very close to each other differing only 
randomly from an ideal shape representative for the 
object kind. Hence, in such cases the mean contour is 
indeed a reliable representative of the undeformed body 
shapes of objects of the same kind. In this sense, the 
classification problem encountered in the present work 
is not a trivial one. All considered "objects" are parasites 
of the same kind and, as Figures |2lJ [22] indicate, their 
undeformed shape contours and the class representative 
mean curves are very close to each other making parasite 
family separation a difficult task. So, in order to obtain 
a high parasite family identification rate, we have de- 
veloped a tuning process for the family representative 
curves in order to increase the separability between the 
corresponding groups. 

8 Conclusion 

In this paper, a new methodology has been introduced 
that tackled two goals: 

1. To exploit images of elastic body deformation in- 
stances, so as to verify assumptions about its mechano- 
elastic properties. The validity of these assumptions al- 
lows for unwrapping the body, i.e. for obtaining the body 
undeformed version from its deformed image. This has 
been achieved here a) by detecting elastic deformation 
invariant curves (neutral line and its cross sections), b) 
by applying to the image of a 2D body deformation the 
inverse deformation process, as it is constructed by a 
sequence of image operations equivalent to solving the 
differential equations governing the body deformation. 

2. To automatically classify the deformed bodies (here 
parasites) on the basis of comparison of the contours 
of their undeformed version. To achieve this, a curve 
classification method has been developed that creates a 
representative curve for each group of the Training Set 
and classifies each element of a Test Set according to 
the distance of the unwrapped contour of this element 
from all these representative curves. Application of this 
methodology to 317 images of highly deformed parasites 
of domestic animals, offered straightened contours and 
body versions that seem to be consistent and reliable 
representations of the undeformed parasites. In addition, 
we have obtained the undeformed shape of fibers, cells 
and human lips by application of the methodology 
introduced here, so as to make its applicability to various 
entities more clear. 

Concerning the problem of deformed parasites identi- 
fication, which was the motivation of the present work, 
we would like to point out that the expert parasitologist 
could classify the gene from the obtained parasite un- 
wrapped versions. However, the authors have employed 
these unwrapped parasite versions developing an auto- 
matic classification system for the parasite genera. The 
classification method developed by the authors, succeeds 
in classifying the deformed parasites to proper groups 
and six families with more than 97.6 % success rate for 
200 randomly chosen combinations of Training and Test 
Sets. 



Appendix A 

Here, we analytically construct the solution of the PDE 
problem 



d_ 
ds 



_9_ 
dn 



(ln||Vw|| 



(31) 



with initial conditions {w(x, y)} (0, S) = M$~ [S (x, y)] 
and 1 0o (^,2/)} (0,5) = (In || Vw\\) (0, S). For solving it, 
first we express the directional derivative of w n at an 
arbitrary direction v, 



v T H(w)h = w n iv T l + w nn v T h 



(32) 



But as it has been mentioned in Sect. 13.31 curvature 
is related with functional w via the relation 0. 
Performing dot product of ^ with h, there results 

™ nl ~™J™ nn ^nFl = 0. By substituting in 

|32| we obtain w nv = w nn v T n. This relation bounds 
w nv and, consequently, Wnv — 



\\Vw\\ 



= £(ln||Vt5||). Writing 



(ln||Vw||) as an extreme of (In ||V^|[ 
obtains: 

For s : k(s) < 0, 

-^(ln||Vt5||) = ^ sup {t^V(ln||Vt5||)} = 

\\v\\<ds 

sup [{ln||Vw||((^, 2 y)+i;)}( S ,5)]-{ln||Vw||(x, 2 y)}( S ,5) 



one 



= lim J! 

ds^O 



OS s 



ds 

& sup [{ln||Vt5||((a;,y)+iO}(0,(S)] 



\\v\\<s 



dilation 
scale of s. 

For s : k(s) > 0, 

-&(ln||Vt5||) = 



[Mx,y)}(0,5)] r i. 
of [Mx,y)}(0,5) 



e. the derivative of the 

{In || ViS|| (*,!/)}(<>,*) at 



Ts^ KV(ln||V^||)} 

\\v \\<ds 



lim 

ds^O 



inf [{ln||Vw||((x, 2 y)+i;)}( S ,(5)]-{ln||Vw||(x, 2 y)}( S ,5) 

\v\\ <ds 

ds 



f s inf [{ln\\Vw\\((x,y)+v)}(0,6)} 

\\v\\<s 

|^| 0o 2/) | (0,5) j, i.e. the derivative of the 
erosion of j^oOz, 2/)} (0, S) = {In || Vw\\ (x, y)} (0, S) at 



as i 



scale of s. 

Then IpT) implies 



d_ 
ds 



M+ 



Inw 



ds 

4-M. 



[Mx,y)j(0,S) 
\Mx,y)\(0,S) 



k(s (5)) < 

k(s (S)) > 

(33) 

where sq(S) is given for r5-deformation using ( [13} and 
the initial condition |0o(^ 5 2/)j (0, S) can be derived 
immediately from a dilated version of the global ini- 
tial conditions (so,#o)- To unify both these cases in 
one process we define the function {n^x, y)} (s (5)) = 
-sgn(«(s ((S)))oo. 
We adopt the process 

a[g(5)](a) = sup [M"[ 5 (5)],inf [M+[g(6)}, (5))]] 

(34) 
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The derivative of this functional reads 
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d_ 



a\g(S)](a) = 



£M-[g{5)} 



^(s (S))<M+[g(S)]^m 
k(s {S)) < 

k+(8 (5)) > M+[g(S)] & 
k(s (S)) > [2] 
(35) 

|0o(x,^)| (0, 5) filter a on s-scale 
space gives a solution for (33} and equivalently of pi) . 
Finally, correspondence between ( [31) with initial condi- 
tions g(S) and filter a[g(S)](s) results {w(x, y)} (s, S) = 



Letting g(S) 
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dru> = llVull — w 



dsl(s) +dsS(s) \\\7w\\h(s) 



(36) 



dr L , = \\Vu\\ +w 



\Vw\ 



dsl(s) -dsS(s) \\Vw\\h(s) 



(37) 



(?u> - \\Vu\\ Z(s)) • D'A' = f^/ • £>'A' = S(s) \\Vw\\ (38) 

(Vz/ - ||Vfi|| Z(s)) = F L / = ||Vw|| (39) 

It also holds that 



d _ 

ds U 



\Vu\\i{s) 



^fu-\\Vu\\i{s) 



(40) 



But || Vu\\ l(s) is the tangent vector of the deformed sym- 
metry line at length s, namely L^(s) = r M (s) = \\Vu\\ l(s). 
Now, if we let £ be the angle between rjj' — L^(s) and 
DM' and 9 the one between — L^(s) and D'i', then, 
combining (381, {39) and (40> we obtain ^'' D '^, = ^rf = 

° r L r-D'A' cost/ 

— 1 and, as a consequence, cos £ = — cos 9, which results 

£ + 9 = 7T. 
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